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Abstract 

The Hamiltonian formulation of the motion of a spinning relativistic particle 
in an external electromagnetic field is considered. The approach is based on the 
introduction of new coordinates and their conjugated momenta to describe the 
spin degrees of freedom together with an appropriate set of constraints in the 
Dirac formulation. For particles with gyromagnetic ratio g = 2, the equations 
of motion do not predict any deviation from the standard Lorentz force, while 
for (7 7^ 2 an additional force, which corresponds to the magnetic dipole force, 
is obtained. 



The description of the motion of a spinning relativistic particle in 3+1 dimensions 
has been a subject of research for many years. Equations describing the evolution 
of the classical spin in a uniform and static external electromagnetic field have been 
proposed by several authors |l|, Q . Modifications of the Lorentz force law for spinning 
particles moving in a slowly varying external electromagnetic field together with some 
generalizations for the equations describing the spin precession in this case have also 
been suggested Some Lagrangian and Hamiltonian formulations, which allow 

one to obtain the equations of motion from an action principle, have been proposed J^- 
T3| . Most of them rely on the use of Grassmann variables |01-[|lOl- In IHl the analogy 



with the relativistic quantum equations was used as a guiding principle. Recently, 
in Ref. []T2[ a canonical procedure has been proposed based on a Routhian function 



Ij] from which the equations of motion for a spinning relativistic particle coupled to 
an external electromagnetic and a gravitational field can be derived. Also, in |jl3l a 
Lagrangian formulation based on classical commuting variables was proposed. It is 
however amusing that there is, even today, an ongoing debate as to the torque and 
force on such particles. 

In this letter we shall consider the Hamiltonian formulation for the description 
of the motion of a spinning relativistic particle in 3+1 space-time dimensions in an 
homogeneous constant electromagnetic field. The model is constructed by analogy 



with the one presented in [|T^] for anyons in 2+1 dimensions. It uses the idea of 
expressing the antisjTumetric spin tensor in terms of some new coordinates and their 
conjugated momenta as in the case of the orbital angular momentum tensor. We start 
with a Hamiltonian formulation of the problem by introducing the corresponding 
constraints in the theory. The Lagrangian can then be derived by using the usual 
Legendre transformation. We hope that the present model will provide an alternative 
formulation to the ones discussed in the literature and mentioned above. 

Let the position of the relativistic particle be described by the four vector x'^ 
(/i = 0,1,2,3), and the spin of the particle by the antisymmetric spin matrix S^^^. 
We shall introduce the auxiliary coordinates n'^ to describe the spin. Therefore, in 
the phase space we work with the configuration space variables x'^ , n^^ and their 
conjugated momenta , p^^\ We can write the following Poisson bracket (PB) 
relations: 

{x^ , p,y} = -g^,u , {Xf^ , x„} = {p^ , p^} = , 
{n, , = -9,. , {n, , n,} = {p^^^ , p^")} = , (1) 

where the metric tensor g^j^^ = diag{l, —1, —1, —1). 
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Let us define the antisymmetric spin matrix Sfj,u as follows: 

S,. = n.p^J^^ - nj;^ . (2) 
The total angular momentum matrix M^^, can thus be written as 

Mfj^i, = Lfj_y + Sfj,ii , (3) 

where L^^ = x^p^ - x^p^ . 

It is straightforward to verify that the four-vector p^ and the tensor M^,^ satisfy 
the Poincare algebra 

{M^^ , Px] = -g^,\Pu + g^xp^, , 
{M^, , Mxf,} = -g^^xM.p - g.pM^x + g^^pM^x + gvxM^p . (4) 

Now we would like to consider the interaction of a spinning particle with gyro- 
magnetic ratio g = 2, electric charge e and mass m with a uniform static external 
electromagnetic field. A vector potential describing such a field can be expressed as 

A, = -\f,.x'' , (5) 

where F^^ is the electromagnetic field tensor. 

In order to describe the motion of such a particle we shall follow a Hamiltonian for- 
mulation analogous to the one presented in [|l^] and based on the imposed constraints 
and the algebraic properties of the system. 

Let us consider the following extended (with respect to the Poincare) algebra: 



{M^y , tta} = -g^j^xT^v + guxT^p , 
{M^^ , Mxp} = -g^xMyp - gypMp^x + gp,pMyx + g^xMpp , 

{Fp. , tta} = , 
{Fpu , Fxp} = , 

{Mpy , Fxp} = —g^xFup — gupFpx + g^pFux + guxFpp , (6) 

where TTp = Pp + eAp. 

It is not difficult to prove that the operator 

C = vr^ - eFp.M^^'' (7) 
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is the Casimir operator of the algebra (|]). By using the relation 

F^,L^^ = 4(M) , 
Eq.(0) can be written in the form: 



with TTfj^ = Pfj_ — eAfj^. By fixing the value of C to be equal to we can then take the 
quantity 

$ = TT^ - eF^.S'"' - (9) 



as a first-class constraint ||T6[ in the Hamiltonian formulation. 

Next we should require that in the rest frame of the particle the spin tensor S^^p has 
only three independent components. This is usually done by imposing the conditions: 

^M-i" = , (10) 

where is the four-velocity and the dot denotes the derivative with respect to the 
proper time r. Since we are considering the Hamiltonian formalism we should employ 
the momenta instead of the velocities. In an external electromagnetic field we thus 
expect the use of vr^ instead of x^. On the other hand, it is easy to see that from the 
definition of the spin tensor (H) it is sufficient to impose the constraints 

(7rp(")) = , (vm) = , 

in order to guarantee that 

S^.vr'^ = . (11) 

Therefore, in our Hamiltonian formulation we shall consider the following set of 
primary constraints: 

$ = TT^ - eF^pS^'' - , 
^1 = (vrp^")) , 

^2 = (™) • (12) 



Since 



{$,9^1} = 0, 

{$,</'2} = 0, 

{^,,^,} = n' + ^F^^S'^'', (13) 
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the constraint $ is a first-class constraint, while and form a pair of second-class 
constraints. We shall define the total Hamiltonian of the system as 

H = + Ai(/?i + \2^2 , (14) 

where A, Ai and A2 are Lagrange multipliers and the canonical Hamiltonian H^an is 
taken to be equal to zero. 

From the consistency conditions on the primary constraints (0), i.e. from the 
conditions 

<i> = {<i>,if} = o, 

</3i = {ypi, H} = A2(7r2 + ''-F^^Sn = , 
= {</^2, H} = -Ai(7r2 + ^-F^^S'^') = , (15) 

it follows that 

Ai = A2 = (16) 

on the constrained surface and that no new constraints appear in the Dirac algorithm. 
The Lagrange multiplier A remains undetermined. 

Notice that the number of physical degrees of freedom in configuration space is 
equal to 8 — 1 — 2/2 = 6. It is thus expected that 3 of these correspond to the position 
of the relativistic particle and the other 3 to the components of the spin. 

It is useful to define the spin four-vector through the relation: 

1 tt" 

= 2^ mP-y—l=^S^"' , (17) 



which in the absence of an external electromagnetic field reduces to the usual Pauli- 
Lubanski four- vector |17]. Using Eqs. (|ll|) and ([l7|) it is not difficult to prove that 
the spin matrix S^y can be written on the constrained surface as follows: 
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It also follows that 



= S^S^' = -^S.^S^-' . (19) 
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The Dirac-Hamilton equations of motion which follow from the PB (|l|) and the 
Hamiltonian (|Tip read 

= 2AeF^,p^^> , (20) 



and they should hold together with the constraints (jT^)- 

Let us now rewrite these equations in a more familiar way. First notice that 

g 

Using Eqs. (|20|) and choosing A = l/2m we obtain 

m 

Also from the definition of the spin matrix (cf. Eq.(0)) and the equations of mo- 
tion (^) it follows that 

S^u = F'^^ {S^aQvP + Sypg^a) . (22) 



Eqs.(pTD and ( p2D imply that the Frenkel conditions (|1^) are preserved in proper 
time, i.e. 

^ {S,,x') = . (23) 
ar 

Also from (^) and (^2]) it follows that the value of the spin is a constant of motion: 

d_ 

dr 



= . (24) 



Finally, Eq.(^) can be rewritten in terms of the spin vector defined in ([T7|). We 
have 

5'^ = -F^uS"" , (25) 
m 

which coincides with the Bargmann-Michel-Telegdi equation for spinning particles 
with g = 2. 

For completeness, let us derive the Lagrange function corresponding to the total 
Hamiltonian (|14D . The derivation is very similar to the one presented in [|l^ . The first 
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step is to express the Lagrange multipliers A, Ai and A2 in terms of the generalized 
coordinates and velocities. This can be achieved by using the relations 

= 2KeF^,n' + Aivr^ , (26) 

and the constraints (|12D. 

The Lagrange function is equal to 

L = (xp) + (np^'^y) = (xtt) + (np^")) + e{Ax) , (27) 

where the momenta tt^ and p^^^ should be expressed in terms of the generalized 
coordinates and velocities using Eqs.([26|). Since 



then we have 



L = 2m^A { 1 + ^F'^uS^" ) + e{Ax) . {21 



Up to terms linear in the electromagnetic field strength we find from Eqs. (|26|) and 
(0): 

A =-^^ + 0(6^^2) , (29) 



and thus we obtain finally 

L = (1 + ^i^M.^'^^) + e{Ax) + 0{e'F') , (30) 



where 



n2 ^ ' (n2)3/2 



It is worth mentioning that in the absence of an external electromagnetic field 
from the equations of motion (pO]) it follows that n'^ = 0. In this case it is more 
convenient to adopt a first order formalism {e.g. as advocated in |T^, [1^) and take 
the Lagrangian as 

L = {xp) + {np^'^y) - A(p2 - _ X^{pp(^)) - X^{pn) . 
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We have presented the Hamiltonian description for the motion of a particle with 
gyromagnetic ratio = 2 in an external electromagnetic field. The main ingredients of 
the model are the following: the phase space of the system is given by the coordinates 

and ra^ and their conjugated momenta and p^"^ respectively; the antisymmetric 
spin matrix S^^j is defined in terms of the four vectors and p^"^ (see Eq.(^) by 
analogy with the orbital angular momentum tensor; the total Hamiltonian is given 
in ([l^) with the primary constraints defined in (|T2|). Out of these starting points we 



have obtained the equations of motion (|2T| ) and (^) (or equivalently Eq. (pS])). 

Let us now consider the case of a particle with gyromagnetic ratio g ^ 2. A 
straightforward generalization of our previous {g = 2) analysis to the case g ^ 2 
would be the introduction of the factor g/2 in the interaction term eF^^^S^^ of the 
constraint $ in Eqs.(p^. However, in order to keep the latter constraint as first-class 
and the constraints yji and (^2 as second-class, we should further modify the constraint 
$. Taking the linear combination 



and imposing the conditions {^,ipi} = {^,ip2} = for the constraint $ to be first- 
class, we can determine the coefficients A and B. Finally we find 

i . _ _ „^ _ _£(|Z^^»..„5,,.^ . (32) 



The equations of motion which follow from the new Hamiltonian 

H = + Ai(^i + A2<^2 (33) 
(with (fi and (p2 defined in (|T^) read as follows: 

™m = 2A71'' ^^^^ ~ F «7r"5^'' 



F^.Sf^'^y 



P^;^ = AegF,.pi-> - ^^,^;^f^J,.^ F-^.^^vr, . (34) 
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For the spin tensor 3/^1,^ and vector we get in turn: 



Ae{g - 2) 



vr 



= AegF^yS" - ^^^^^^tt^F.^vt"^^ 



We note that from the equations of motion (^)-(^) it follows that the value of 
the spin (|19D is preserved in time. Moreover, using the expression for tt'^ which follows 
from the first equation in (|34D , i.e. 

2A\' ^ {rr' + IF.yS^'^r''' ) ' 
and substituting it into the definition of the spin vector (IT^), it is easy to show that 

S, = —^e.^p.x'^S^^ . (36) 



Thus the condition S'^x'^ = is automatically satisfied by the equations of motion. 

For the case 5^ = 2, if one chooses A = l/2m, Eqs.(p4D-(p5D coincide with Eqs.(pl|) 
and (p5l). However, when g ^ 2, additional quadratic and higher order terms in 
the field strength F^^, appear. These terms reproduce the Bargmann-Michel-Telegdi 
equation 0] in the linear approximation. Indeed, from (|3^ ) - (|35|) we obtain 



X, = -F,,x^ + {g- 2)0{e'F^) 
m 



S, = ^F,yS^ - F^ ^^S^ + {g- 2)0{e'F') . (37) 

2m 2m 



We remark that the extra terms of order e^F^ in Eqs. (|37|) do not coincide with 
the ones previously obtained in the literature 0, 0, |13|. One way to see that is 
to consider the nonrelativistic limit. For particles with g = 2 our equations do not 
predict any deviation from the standard Lorentz force, while in Refs. [^, [T^ a 
term proportional to gB x (S x E) is obtained. For g ^ 2, however, Eqs.(|3^) do 
predict in the nonrelativistic limit an additional force, namely, 

-%^-'-^)' 

which modifies the Lorentz force. 

We have considered so far the case when the electromagnetic field strength F^^, 
is constant. We can modify the model to obtain the equations of motion in slowly 
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varying fields. In the latter case, the constraint $ defined in Eq.(^) as well as $ in 
Eq. ( P^D will also generate the correct equations of motion so that in the nonrelativistic 
hmit we obtain the term ^V{B ■ S) §, |T2|, |13 . 



We are grateful to M. Bander, I. Khriplovich and G. Kunstatter for useful discus- 
sions. 
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